Neutral Color Superconductivity Including Inhomogeneous Phases 

at Finite Temperature 
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We investigate neutral quark matter with homogeneous and inhomogeneous color condensates at 
finite temperature in the frame of an extended NJL model. By calculating the Meissner masses 
squared and gap susceptibility, the uniform color superconductor is stable only in a temperature 
window close to the critical temperature and becomes unstable against LOFF phase, mixed phase 
and gluonic phase at low temperatures. The introduction of the inhomogeneous phases leads to 
disappearance of the strange intermediate temperature 2SC/g2SC and changes the phase diagram 
of neutral dense quark matter significantly. 
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I. INTRODUCTION 

The emergence of gapless color superconductivity [l|- 
4] promoted great interest in the study of dense quark 
matter 0, Q. A crucial problem is whether the two 
flavor gapless (g2SC) or breached pairing[7i] color su- 
perconductivity is stable. The first instability is the 
thermodynamical instability when the chemical poten- 
tial difference between the two species is fixed. This is 
called Sarma instability Q. It is now widely accepted 
that the Sarma instability can be cured under charge 
neutrality constraint due to the long range electromag- 
netic gauge interaction [H Q. The second instability 
is the chromomagnetic instability described by nega- 
tive Meissner mass squared [§-[l3[ or negative superfluid 
density [Til. Il5j. In this case, some inhomogeneous con- 
densed phases are energetically favored [16|-|23|. A specific 
character of the chromomagnetic instability is the nega- 
tive Meissner masses squared for the 4- 7th gluons in the 
gapped 2SC phase [24j. The third instability is the insta- 
bility against the phase separation or mixed phase when 
the surface energy and long range gauge interaction are 
excluded. 

All these instabilities are confirmed only at zero tem- 
perature and in the weak coupling region. While it 
is recently argued that the breached pairing phase will 
be free from negative superfluid density or Meissner 
mass squared in the strong coupling BEC region(25l- 
[28j . the diquark condensate at moderate baryon den- 
sity seems not in this region. As for the tempera- 
ture effect, in the study of isospin asymmetric nuclear 
matter |29|, |3(|, two flavor color-superconducting quark 
matter[3| and general breached pairing superfluidity [3l| . 
the fermion pairing correlation behaviors strangely: 
The superfluid/supercondutor order parameter is not a 



monotonously decreasing function of temperature and its 
maximum is located at finite temperature. Especially, for 
large enough density asymmetry between the two species 
at fixed coupling or for weak enough coupling at fixed 
density asymmetry, the superfluidity appears at inter- 
mediate temperature but disappears at low temperature. 
This strange temperature behavior is quite universal and 
also discussed in the study of atomic Fermi gas with den- 
sity imbalance [32H341 ] . It is found [H, HH that the super- 
fluid density and number susceptibility of the uniform su- 
perfluid phase are positive only in a temperature window 
near the critical temperature. From the recent instability 
analysis at finite temperature 34] , if LOFF phase is taken 
into account, this strange temperature behavior will be 
washed out and the superfluid order parameter remains 
a monotonous function of temperature. Therefore, for 
the color superconductivity at moderate density, we may 
have the following estimations: (1) The superconductor is 
free from the chromomagnetic instability and also stable 
against the mixed phase at temperatures near the criti- 
cal value. (2) The phase diagram of neutral dense quark 
matter is significantly changed if some non- uniform phase 
such as LOFF is taken into account. In this paper we will 
show that these two estimations are true. 

The paper is organized as follows. In Section [TTJ we 
study a naive two-species model which possesses the basic 
mechanism for positive Meissner mass squared and gap 
susceptibility at finite temperature. We then calculate 
in an extended NJL model the Meissner masses squared 
for gluons and the gap susceptibility in two flavor color- 
superconducting quark matter with charge neutrality in 
Section IIIII and investigate the LOFF phase and gluonic 
phase with charge neutrality in Section IIVI We summa- 
rize in Section IVl 



II. RESULTS IN A TOY MODEL 
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In this section we study a toy model containing two 
species of fermions proposed inflj, [H, [2(|. While this 
model only reveals the chromomagnetic instability for the 
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8th gluon and the important non-abelian type chromo- 
magnetic instability is excluded, it can help us to better 
understand why the choromomagnetic instability can be 
cured at finite temperature. The Lagrangian of the model 
is defined as 

£ = ^(i7^-m)V + G(^ c i75TiV) (#75^°) , (1) 

where ip = (ipii'4 , 2) T and ip — (f/'i,^) are Dirac spinors 
including two species of fermions, tp c = C^p T and rp — 
ip T C are charge-conjugate spinors, C — ij 2 j° is the 
charge conjugation matrix, m is the fermion mass, G 
is the attractive coupling, and T\ is the first Pauli ma- 
trix in the two-species space. In the following we take 
m = which corresponds to the quark matter at high 
baryon density. It is quite convenient to work with the 
Nambu-Gorkov spinors defined as <I> — (ip,ip c ) T and 
^ = ip c ). We assume that the attractive coupling 
is not very strong and the BCS mean field theory re- 
mains a good treatment at finite temperature. The study 
on possible BCS-BEC crossover in relativistic fermion su- 
perfluid can be seen in (35l . |36| . In the mean field approx- 
imation, the pair fields are replaced by their expectation 
values which serve as order parameters of superfluidity. 
Introducing the order parameters 

A = -2G (^ c i 15 t^) , A* = -2G {^i l5 T^ c ) (2) 

and taking them to be real, we can express the thermo- 
dynamic potential as 



4G 2 ^ 



d 3 p 

(27T) 5 



TrlnS 1 (uj n ,p) (3) 



with the mean field fermion propagator in the 4- 
dimensional Nambu-GorkovC^two-species space, 



[G+]- 1 nsAn 
i l5 A Tl [Go"]- 1 



(4) 



where [Gg 



(iio n ± /i)7o — 7 • p are the inverse of 



free fermion propagators with frequency uj n — (2n + 
l)irT, n = 1,2,.... We have introduced the chemi- 
cal potentials /ii and \i 2 for the two species via adding 
the term = Tpftj ^ to the original Lagrangian with 
fx = diag(p\, fi 2 )- After a straightforward calculation, 
the thermodynamic potential can be evaluated as 



where we have defined the notation E^ = 
yj{p± p) 2 + A 2 and introduced the average chemi- 
cal potential £l = + /i2)/2 and the chemical potential 
mismatch S/i = (/i2 — /i2)/2. The gap equation for A can 
be derived via minimizing the thermodynamic potential, 



dn/dA = 



(7) 



and the fermion numbers are determined by the defini- 
tion, 



ii = —dtt/dni, n 2 — —dfl/dfi 2 - 



(8) 



The order parameter A under constraint of fixed numbers 
ni and n 2 can be calculated by solving the gap equation 
§7$ together with the number equations ((H|). Since the 
gap equation for A suffers ultraviolet divergence, we in- 
troduce a three momentum cutoff A to regularize it. 

To model the two flavor dense quark matter with beta 
equilibrium and charge neutrality, we choose the cou- 
pling constant G < 0.1 Go with Go = Air 2 / A 2 to ensure 
weak coupling and fix the averaged chemical potential 
/2 = 0.615A and the number ratio n\/n2 = 2. In Figfl] 
we show the temperature behavior of A and S/i. For the 
coupling G = 0.095Go, the gapless state satisfying the 
condition <5/i > A appears at low and high temperature 
regions but disappear at intermediate temperature. For 
weaker couplings G = O.O8G0 and 0.073Go, the super- 
fluid is always gapless. For weak enough coupling such 
as G = 0.073Go, the superfluidity appears only at inter- 
mediate temperature. 

For a symmetric system with n± = n 2 , it is well-known 
that the temperature effect destroys the pairing, and 
the pairing gap is a monotonously decreasing function of 
temperature. However, for an asymmetric system with 
ni 7^ n 2 j the two species of fermions have mismatched 
Fermi surfaces, and the temperature effect not only de- 
forms and reduces the Fermi surfaces which melts the 
gap, but also makes the overlap region of the two species 
wider which favors the pairing. The competition of the 
two opposite effects results in a non-monotonous temper- 
ature behavior of the pairing gap. This is the reason why 
in the low temperature region in Figfl] the gap increases 
as the temperature increases and even disappears when 
the coupling is small enough. 



Meissner Mass 



A 2 * 

4G ^ 



1=1 



d 3 p 



W{E!) 



(5) 



with W(Ej) = Ei + 2Tln (l + er EljT ), where the sum 
runs over all fermionic quasiparticles. The quasiparticle 
dispersions Ej(p) calculated via det S^ 1 (p) = are given 
by 
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Let us assume that the fermions couple to a U(l) gauge 
field with coupling g, the superfluid becomes a su- 
perconductor where the gauge field obtain the so called 
Meissner mass. The Meissner mass can be calculated 
from the polarization tensor 



IP" (A) 



d 3 p 
(2^)3 



Tr 



f»S(p)r'S(p-k) (9) 



with the vertex r M = 177^0-3 (8 To, where 03 is the third 
Pauli matrix in the two-species space, tq is the identity 
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FIG. 1: The pairing gap A and chemical potential mismatch 
S/i as functions of temperature T for three coupling values. 
A, 8/1 and T are all scaled by the cutoff A. 



matrix in the Nambu-Gorkov space, and the propaga- 
tor S can be obtained from its inverse (@| by taking the 
technics used in [l5j]. The Meissner mass tua is defined 

as 



i\ = - iim(<% - fc, vi k l 2 )n lJ (fc = 0, k). (10) 

Z p— >0 



The calculation of the Meissner mass is quite straight- 
forward but tedious, we quote here only the final result, 



m\ = m 2 , + m 2 , 



(11) 



the diamagnetic term m 2 and paramagnetic term m 2 are 



explicitly expressed as 
2 4 9 2 f d 3 p 



(2tt) 3 



X -(p)(f(E 1 )+f(E 2 )-l) 
X + (P) (f(E 3 ) + f(E 4 )-l)] 



d 3 p 



52 I'm, 



m 2 A 
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where f(x) is the Fermi-Dirac distribution function, f'(x) 
is its first order derivative, and the functions x ± (p) are 
defined as 



X ± (P)=± 



E^Ev 



Et {(EX? 



2pE 



(13) 



with E = p ± fi. Since the diamagnetic term suffers 
ultraviolet divergence, it contains a subtraction term 



m\ 



3 



d 3 p 1 g 2 A 2 



(2tt) 3 p 3tt 2 



(14) 



to ensure m\(A = 0) = 0. 

Above the critical temperature T c where the pairing 
gap A vanishes, we have x (p) = l/(2f>) and 



3 
9 2 



d 3 p 1 ^ 



{^Ypfi 

d 3 p 



f(p-Vi) + f(p + Mi) 
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Taking partial integration we find 



d 3 P fl( , , 2 



d3 P 1 f ( 4- ^ 

;-/(P±M) 



(2tt) 3 p 



(15) 



(16) 



for any T and /z. Therefore, the total Meissner mass m\ 
is zero in the normal phase, as it is expected. 

In weak coupling with A <C M an d ^M ^ Mj the Meiss- 
ner mass squared at zero temperature can be well ap- 
proximated asQ 



5V 
3tt 2 



1 - 



5//6>(5m - A) 
V V - A 2 



(17) 



In the gapless phase with A < S[i the negative Meissner 
mass squared is quite different from m 2 A = g 2 p, 2 /3n 2 > 
in the symmetric BCS case where there is no paramag- 
netic contribution. 

Now the question is: Is the Meissner mass squared 
always negative below the critical temperature T c ? To 
answer this question, we derive the critical behavior of 
the Meissner mass squared at T c . Since the gap equation 
for A below T c contains only A 2 and A continuously 
approaches to zero at T c , we have the critical behavior 
for the order parameter 



A(T) cx (1 - T/T c 



,1/2 



T^T C . 



(18) 
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As a consequence, the Meissner mass squared behaviors 
as 



i\{T) oc 1 - T/T c > 0, T -> T c 



(19) 



Combining m 2 A < at zero temperature and m 2 A > 
for T approaching to T c , m A should be negative at low 
temperature but positive in a temperature window close 
to T c . 

In Figj2] we show the temperature behavior of the 
Meissner mass squared for three coupling values. In any 
case, there exists a temperature window close to the crit- 
ical temperature where the Meissner mass squared m A is 
positive which means that the gapless phase is magneti- 
cally stable and hence stable against the LOFF phase. 
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B. Density Fluctuations 

The general stability condition against changes in num- 
ber densities for a two-component system is described 
by the total free energy of the svstemfdH [421 ], F = 
J d 3 xJ r [n cr (x)](cr = 1,2). Considering its fluctuations in- 
duced by small number changes 8n a (-x.), the first-order 
variation 6F vanishes automatically due to the charge 
conservation, J d 3 xSn a (x) = 0, and the second-order 
variation S 2 F is expressed in the quadratic form 



6 2 F = 



1 



dh 



E 



d 2 T 
dn„dn c 



-5n a Sn 



(20) 



Therefore, to achieve a stable homogeneous phase, the 
2x2 matrix d 2 J r /dn a dn a ' should be positively definite, 
namely, it has only positive eigen values. From the rela- 
tion between the free energy J-(n a ) and thermodynamic 
potential £l(/i CT ), Q = J- + /Jini + /i2«2, it is easy to 
check that the stability condition to have positively def- 
inite matrix d 2 T /dn a dn a ' is equivalent to the condition 
to have positively definite number susceptibility matrix 
— d 2 ^ljdpL a d[i a i . 

For systems without mass difference between the 
two species, the condition to have positive eigen- 
values of —d 2 Q/d[i a d[i r7 i can be reduced to the 
condition [25j that the imbalance number susceptibility 
X = —(d 2 n/dSfi 2 ) tl = (dSn/dSn)^ should be positive. 
For x < 0, the density difference Sn — ri\ — ni increases 
with decreasing chemical potential difference Su, which 
is certainly unphysical and means that the uniform phase 
is unstable against density fluctuations. Employing the 
gap equation which determines the condensate as a func- 
tion of chemical potentials, we can express the imbalance 
number susceptibility x as a direct and an indirect part, 



_ fdSn\ /<95n\ / dA 



(21) 



From the expression 
d 3 p 



Sn 



- f(E 2 ) + f(E 3 ) - f(Ei)] (22) 
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FIG. 2: The Meissner mass squared m\ scaled by m% = 
g 2 fl 2 /3n 2 and the gap susceptibility ka scaled by kq = 2fi 2 /it 2 
as functions of T/A for three coupling values. 



Since (dSn/ 'dSu)^^ is always positive, the stability con- 
dition x > is controlled by the gap susceptibility 
ka = (d 2 fl/dA 2 )p / which determines if the solution of 
the gap equation is the minimum of the thermodynamic 
potential. 

The gap susceptibility ka can be explicitly evaluated 
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ka = 2 



d 3 p 
(2tt) 3 



A 2 /l-/(£a)-/(£ 2 ) 



A / \ ^A 

+f'(E 1 ) + f(E 2 ) 
A 2 /l-/(£ 3 )-/(£ 



£1 



AJ \ 

+/'(^ 3 ) + /'(£*) 



(24) 



where we have considered the gap equation for the con- 
densate A. At weak coupling and at T — 0, ka can be 
evaluated as 



2fl 2 



Sfid{Sn - A) 
- A 2 



(25) 



which leads to x < in the gapless phase. Note that 
while the Sarma instability ka < can be cured via 
charge neutrality constraint, the instability ka < in- 
duced by density fluctuations can not be removed by 
charge neutrality. Similar to the Meissner mass squared, 
ka near T c takes the form ka oc 1 — T/T c > 0. There- 
fore, the uniform superfluid phase is stable against den- 
sity fluctuations at temperature close to T c . 

The temperature behavior of the gap susceptibility ka 
is illustrated in Figf21 In any case, there exists a temper- 
ature window close to the critical temperature where ka 
is positive which means that the gapless phase is mag- 
netically stable and hence stable against the phase sep- 
aration. Note that the stable region against the LOFF 
phase is larger than the stable region against the phase 
separation. 



III. NEUTRAL 2SC/G2SC PHASE 

We investigate in this section the two flavor color su- 
perconductivity in an extended NJL model. The La- 
grangian density of the model including quark-quark in- 
teraction sector is defined as 



C 



ip (ij^dfj. - mo) V + G s 



+G D (^i^ee^ip) (ipi^ee^ip ) , 



where Gs and Go are, respectively, coupling constants 
in color singlet and anti-triplet channels, the quark field 
ipia with flavor index i and color index a is a flavor dou- 
blet and a color triplet as well as a four-component Dirac 
spinor, t = (ti,T2,T3) are Pauli matrices in flavor space, 
and (e)ij = e iJ and (e 7 ) Q /3 = e a ^ 7 are, respectively, to- 
tally antisymmetric tensors in flavor and color spaces. 
We focus in the following on the color symmetry break- 
ing phase with nonzero diquark condensates defined as 

7 - -2G D (^ c i^se^}, 

0* = -2G c (^7 5 £e^ C ), 7=1,2,3. (27) 



To ensure color and electric neutralities, one should intro- 
duce a set of color chemical potentials fi a {a — 1, 2, 8) 
with respect to color charges Qi, Q2, Q% and an elec- 
tric chemical potential /i e with respect to the electric 
charge Q e . The ground state of the system is deter- 
mined by minimizing the thermodynamical potential, 
dQ/dfiy — under the charge neutrality constraint 
Q e = and Q a = 0(a = 1, 2, 8). 

Since the model Lagrangian is invariant under the 
color SU(3) transformation, we can choose a specific 
color symmetry breaking direction. The most convenient 
choice is 



Si = fa = 



= A ^ 



(28) 



To simplify the calculation, we consider the chiral limit 
with too = and assume chiral symmetry restoration 
in the color superconducting phase. This assumption 
is confirmed when the coupling constant Go in the di- 
quark channel is not too large [33, HI]. In the specific case 
(j28|) . Qi, ■ ■ ■ ,Qi vanish automatically and only Qg can 
be nonzero. Therefore, we can introduce the color chem- 
ical potential fj,s only, and the quark chemical potential 
matrix elements can be expressed as 



a/3 



(flSij - HeQij)$a8 + —7=^sSij(T S ) a p. (29) 
V3 



In mean field approximation the thermodynamical po- 
tential of the system can be expressed as 



CI 
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4G 



T 



D 



E 



d 3 p 
(2^f 



TrlnS" 



(30) 



where £l e is the contribution from the free electron gas 



12tT 2 



Me + 27r 2 M 2 T 2 
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-T 



(31) 



The inverse of the quark propagator S 1 in the Nambu- 
Gorkov space can be written as 



-i7 5 ee 3 A 



[Go]- 1 



(32) 



with [Gq] 1 = (iu>„ ± pff)jo After a straightfor- 



(26) ward algebra Q can be evaluated as 



n = 



/A 2 
AG, 



-5> 



d 3 p 
(27p 



W{E! 



(33) 



where the sum runs over all quasiparticles. The quasi- 
particle dispersions Ej(p) calculated by det S^ 1 (p) = 
are given by (J6j> and 



E 5 (p) 

Mp) 



: P- 
: P- 



E e (p) 
E 8 (p) 



P 
P 



(34) 



with the chemical potential mismatch Sfi = fJ, e /2 and 
averaged chemical potential /2 = fi — fi e /6 + /is/3 for 
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paired quarks, where \i = fJ.s/3 is related to the baryon 
chemical potential hb- The degenerate factor Nj is 2 
for I = 1,2,3,4 and 1 for I = 5,6,7,8. Minimizing the 
thermodynamic potential Q, 



dil/dA = 0, 



and considering the charge neutrality conditions 



dn/dfis = 0, dn/dn e = 0, 



(35) 



(36) 



we can determine simultaneously the order parameter A 
and chemical potentials /i e and fig in the neutral uniform 
color superconductor. At weak coupling, the color chemi- 
cal potential fi$ is only a few MeV, and the electric charge 
neutrality plays the role of the condition m/n^ = 2 in 
the toy two-species model. 

The numerical solutions of the condensate A and 
chemical potential mismatch Sfi are demonstrated in 
Fig 0] at fixed quark chemical potential /i = 400 MeV for 
several values of coupling G n . There are three parame- 
ters in the model. The momentum cutoff A and coupling 
Gs can be fixed by fitting the pion decay constant and 
chiral condensate in the vacuum and we denote the 
coupling Gd by the ratio rj = Gd/Gs- For weak coupling 
such as r\ = 0.75, 0.70, 0.66, the temperature behavior of 
A and Sfj, is similar to that in the toy two-species model. 
The uniform color superconducting phase disappears at 
the critical coupling rj c ~ 0.63. For strong coupling such 
as rj — 0.85, the quark matter is in gapped phase in a 
wide temperature region and in gapless phase only in a 
small window close to T c . Again, the strange behavior 
of the gap in the low temperature region for small r\ is 
due to the competition of the two opposite temperature 
effects for pairings with mismatched number densities. 



A. Meissner Mass 

The Meissner masses for gluons and photon can be 
evaluated via the polarization tensor [9j 



d 3 p 
(2^3 



Tr 



T^ a S(p)nS(p-k) , (37) 



where the vertex in the color-flavor space is defined as 
f £ = diag (g.7"T 8 , -5,7"^) for a = 1, 8 and f £ = 
diag{e^Q,-e^Q) for a = 7. In the 2SC/g2SC phase 
with the conventional choice (|28p. the diquark conden- 
sate breaks the gauge symmetry group SU(3) C ® U(l) em 
down to the SU(2) C ® U(l) em subgroup. Therefore, we 
need to calculate the Meissner masses for the 4-8th glu- 
ons and the photon only. Most of the analytic calculation 
is presented in [9|. With the function A, A', B, C, D, H, J 




10 20 30 40 50 60 
T(MeV) 




10 15 20 25 
T(MeV) 



FIG. 3: The pairing gap A and chemical potential mismatch 
Sfj, as functions of temperature T at fj, = 400 MeV for four 
coupling values in the diquark channel. 
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listed in Appendix lAl we obtain 

/ dpp 2 [A + A' - 2B + 4{C - D)] , (38) 

97T J 

— f dpp 2 UA + A' + 4B + 2(5C + AD)} , 
/■ dRp2 [2A _ A , _ B + 2{c _ D)] 

in J 



with a s = g 2 /(4ir) and a e — e 2 /(4ir) for the 8th gluon 
and photon and 



m A = m AA = m 55 = m Kfi = m 



'66 



''77 



=^y ^ 2 (i/+ 2 j) 

(39) 



for the 4- 7th gluons. Due to the nonzero m| 7 , the 8th 
gluon and the photon mix with each other, and the physi- 
cal Meissner masses squared are given by the eigen values 



2 1 

TO 88 = 2 



'88 
l 88 



77 
,2 



(in 
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.(40) 



When the coupling is not very large, we have approx- 
imately m 2 - ~ 0, which is consistent with the analy- 
sis of symmetry breaking in weak coupling limit. Since 
a s r~j 1 >> a e ~ 1/137, we have 



m h 



'88 ■ 



(41) 



With the explicit form of the functions A,A',B,C,D 
shown in Appendix El m| takes the same expression as 
in the toy two-species model, 



'A ■ 



(42) 



and therefore, vn\ and have exactly the same tem- 
perature behavior. Now the most important task is to 
investigate whether m 2 can be positive at finite temper- 
ature. For T — > T c , we have ^ 8 — > and the functions H 
and J are reduced to 

d 

H = Y,[f(p-IH) + f(p + IH)], J = H/p (43) 



with fj, u = fi — 2fj, e /3 and fid = A* + A*e/3- Employing the 
same trick used in Section HH we find m| = for T > T c 
and 



ml(T) cx 1 - T/T c , T -> T c 



(44) 



for T below but close to T c . Therefore, there must ex- 
ist a temperature window near the critical temperature 
where the two flavor color superconductor is free from 
chromomagnetic instability. 

In FigU]we illustrate the temperature behavior of the 
Meissner masses squared m 2 and m| at fixed chemical 
potential fi = 400 MeV for several coupling values. For 
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FIG. 4: The Meissner masses squared m| and m| scaled by 
m 2 g — Aotsp 2 /3ir as functions of temperature T at fi = 400 
MeV for four coupling values. 



weak coupling 0.62 < 77 < 0.8, there exist two intermedi- 
ate temperatures T4 and Tg , m 2 is negative at < T < T4 
and positive at T4 < T < T c , and m\ is negative at 
< T < T$ and positive at T 8 < T < T c . In a wide 
range of coupling, T4 is larger than Tg. Only for suffi- 
ciently small coupling rj < 0.66, T4 coincides with T 8 or 
even becomes less than T%. For strong coupling rj > 0.8, 
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the 8th gluon is free from chromomagnetic instability at 
any temperature, but the 4- 7th gluons suffer negative 
Meissner mass squared in the low temperature region 
< T < T4 which indicates the instability against the 
gluonic phase [2lj|. 



B. Density Fluctuation 




Like the toy two-species model, we can study the sta- 
bility of color superconductivity against the number fluc- 
tuations of the paired quarks 49] . Similarly, the stability 
condition can be reduced to \ = (9Sn/dSfi)^ > where 
5n = ndi — n U 2 = n c i2 — n u \ is the density imbalance be- 
tween the paired quarks. Using the techniques in Section 
II, the gap susceptibility ka which controls the stability 
takes almost the same form (|24l) , the only difference is the 
replacement of the factor of 2 in front of the momentum 
integration in ([24]) by the factor of 4. 

In FigJS]we show the temperature behavior of the gap 
susceptibility ka for four coupling values in the diquark 
channel. In any case, there indeed exists a temper- 
ature window where ka is positive and therefore the 
gaplcss phase is magnetically stable and hence stable 
against the phase separation. Recently, it is argued 
that the Higgs instability [3^] indicates the instability of 
2SC/g2SC against the mixed phase. The gap susceptibil- 
ity is the long wavelength limit of the Higgs instability. 
For a complete study we may need to check the Higgs 
instability at finite temperature. 



IV. INHOMOGENEOUS PHASES 

We have calculated the Meissner masses squared and 
the gap susceptibility in the neutral 2SC/g2SC phase, 
and found that there exists a temperature window close 
to the critical temperature T c where they are both pos- 
itive. Therefore, the uniform phase is free from the 
chromomagnetic instability and stable against the mixed 
phase in this temperature region. In general case, the 
stable region for the 4- 7th gluons is smaller than the sta- 
ble region for the 8th gluon, namely T4 > Tg. Only for 
sufficiently small coupling 7/ < 0.66, the 4- 7th gluons can 
have larger stable region than the 8-th gluon, namely 
T4 < Tg. For large enough coupling, the 8th gluon is free 
from the chromomagnetic instability at any temperature 
and there exists only a turning temperature T4 where m| 
changes sign. 

If T4 > Tg, the single-plane wave LOFF phase can 
not completely cure the chromomagnetic instability, since 
there exists a region T4 < T < Tg where only m| is 
negative. In this section we will focus on possible inho- 
mogeneous phases in two limits, very weak coupling and 
very strong coupling. For very weak coupling, we have 
T4 < Tg, it has been shown in [22[ that the neutral LOFF 
phase is free from the chromomagnetic instability at least 
at T = 0. For very strong coupling, only 4- 7th gluons suf- 
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T(MeV) 




10 15 20 25 30 
T(MeV) 



FIG. 5: The gap susceptibility ka scaled by kq = 2fj? /n 2 as 
a function of temperature T at /1 = 400MeV for four different 
couplings. 



fers instability, the introduction of gluonic phase may be 
sufficient to cure the instability. 
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A. LOFF Phase 

We firstly discuss the LOFF state. The numerical cal- 
culation will be performed at very weak coupling rj < 
0.66. The order parameters <p 1 and 4>* are complex con- 
jugate to each other and can be set to be real only in a 
uniform color superconductor. In the LOFF phase the 
phase factor of the order parameter is non-uniform. For 
the sake of simplicity we take the following single wave 
LOFF ansatz 



Ae 2i,.x ) 



Ae 



-2iq-x 



(45) 



It has been shown that in the two flavor color supercon- 
ductor the LOFF momentum q can be regarded as the 
8th gluon condensate [22j or a spontaneously generated 
Nambu-Goldstone current [l9j. namely 



1 



With a transformation for the quark fields, 

x{x) = i>(sr)e-«*, x c (f) = v c (zK^, 

we can express the thermodynamic potential as 



n = 



_A2 
AG, 



-E 

2 ^ 



-A^TrlnS, x (iu n ,p) 



(46) 



(47) 



(48) 



with the inverse of the mean field propagator in LOFF 
state 



S 1 {iuj n ,p) 



[G+] 



-i l5 ee 3 A [G-y 



(49) 



where the q dependent free propagators are defined as 



[G, 



±i-i 



(50) 



The neutral LOFF state should be determined self- 
consistcntly by the gap equations for the condensate and 
pair momentum and the charge neutrality constraints, 



an 

OA 



0. 



0, 



0. 



0. 



(51) 



where we have chosen a suitable frame with q — (0, 0, q). 
Since fig, is very small in weak coupling we will simply 
set fis = in numerical calculations. 

Following the treatment in (l6| . we expand the ther- 
modynamic potential in powers of the pair momentum q 
in the vicinity of q = 0, 



9=0 



id 2 n 



(52) 



9 =0 



The linear term vanishes automatically because of the 
gap equation for q. Using the relation 



S- 1 ® =S- 1 (q = Q)-Ttf- 



(53) 



and taking the derivative expansion of in powers of 
the gauge field, we can easily obtain the relation between 
the momentum susceptibility n q and the Meissner mass 
squared for the 8th gluon 



= d 2 n/dq 2 \ q=0 = l2ml/g 



(54) 



On the other hand, since q = must be a solution of 
the gap equation dO,/dq = which corresponds to the 
homogeneous 2SC/g2SC phase, the first order derivative 
of with respect to q must take the form [34[ 



dfl/dq = qQ{q). 



(55) 



The momentum solution for the LOFF state is given by 
Q{q) = 0. The case here is similar to the gap equation for 
the pairing gap A which contains a trivial solution A = 
corresponding to the normal phase and a finite solution 
A / corresponding to the superfluid phase. For the 
formal proof here we do not need the explicit function 
Q{q). From the identity 



8 2 n/dq 2 = Q(q) + qQ'{q) 



and (T54")) . we obtain 



ml = ,g 2 Q(0)/12, 



(56) 



(57) 



which means Q(0) — at T — T§ where m§ changes 
sign. Therefore, the LOFF momentum q must vanish 
at T = Tg, providing that the neutral LOFF solution is 
unique for a given T and y,. This indicates that there is 
no neutral LOFF solution at T > Tg- Below but near the 
temperature Tg, the LOFF momentum q is very small, 
the small q expansion is valid 



n(q) - 0(0) = 6m 2 8 q 2 /g 2 



(58) 



For negative m| below Tg, the neutral LOFF state has 
lower free energy than the uniform state at least at 
T < Tg. Since the small q expansion is like a Ginzburg- 
Landau expansion, we conclude that the LOFF momen- 
tum near the critical point behaviors as 



q(T) ~ (1 - T/T 8 ) 



1/2 



T^Tg 



(59) 



Now we turn to the numerical calculation of neutral 
LOFF state. The explicit form of f2 can be obtained if 
we neglect the mixing between the particles and anti- 
particlesjl^] since we have A, q <§C /2 at weak coupling. It 
reads 



AG, 



J v ; U=l 7=9 



(60) 



where the sum runs over all quasiparticles. The quasi- 
particle dispersions Ej(p) calculated by det^" 1 = are 
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given by 



E 7 ,8(P) 

Eii.uip) 
Ei3,u(p) 
Ei5,ie(p) 



= E7+5E q ±5(i, 



J Aq 
^q 



E^,. - 5E q ± <5/x, 



= E 
= E 



Aq 



Aq ~ 

+ 
W- v\ 
\p + q\ 
\p- o\ 



5 Eg ± <5/i, 

8E q ± (5/1, 

Ms ± Sfi, 
Ms ± <5m, 

Ms ± ^M: 

Ms ± <^M- 



M ' 
A' 
M 
M 



(61) 



where -E AlJ and SE q are defined as 



Aq 



= Tj(\p + q\ + \p-q]-m V4 + A 2 , 



(62) 



Notice that in our ultra-relativistic system, the thermo- 
dynamic potential of the LOFF state suffers from a un- 
physical term ~ g 2 A 2 [22|]. In fact, 0(A, q) should recover 
the result for the free quark gas when A — > 0. Thus to 
study the neutral LOFF state, we define the following 
subtracted thermodynamic potential 



a ub (A, q ) = n(A, q ) - n(o, q ) + n(o, o). 



(63) 



Thus we can set the momentum q to be zero in the dis- 
persion of quasiparticles 9—16 and treat the UV diver- 
gence for the quasiparticles 1 — 8. It is obvious that we 
recover the thermodynamic potential of the 2SC phase 
when q = and the thermodynamic potential of the free 
quark gas when A = 0. We can further make approxi- 
mation on E^ and 8E q since A,g</iat weak coupling 
and all integrals are dominated over the Fermi surface 
\p\ = \x. In this approximation we have 



E 



Aq 
5E a 



Ei 



qcost 



cos 



(64) 



with 8 being the angle between p and q. In this case the 
subtraction term at T = can be analytically evaluated 

as 



Oa = 



d 3 p 
(2+p 



i 



2„2 



2A 2 q 
"3^ 



(65) 



which is from the kinetic term of the quasiparticle disper- 
sion and is indeed diverges as ~ A 2 q 2 . At finite temper- 
ature, however, we do not find an analytical expression 
for J7a- 

In the numerical calculation of the neutral LOFF state 
at finite temperature, we then use the subtracted ther- 
modynamic potential f2 su b(A,q). The gap parameter A 
and pair momentum q as well as the electron chemical 
potential /x e in the neutral LOFF state are determined 



by self-consistently solving the coupled set of equations 
([51]) but with the subtracted thermodynamic potential 
f2 su b- For simplicity, we will neglect the color neutrality 
condition and set p% = 0. 

In Figj6] we display the numerical result of neutral 
LOFF state for -q = 0.66 and 0.70. The neutral LOFF so- 
lution exists only at low temperature, namely at T < Ts, 
and the LOFF momentum approaches to zero contin- 
uously at T — T 8 . Since the LOFF phase has lower 
thermodynamic potential than the uniform phase (we 
checked this numerically), the LOFF phase is energet- 
ically more favored than the uniform superconductivity. 
Considering both the stable LOFF state at low temper- 
ature and the stable uniform superconductivity at high 
temperature, the strange intermediate temperature su- 
perconductivity disappears, and the order parameter be- 
comes a monotonously decreasing function of tempera- 
ture, like in the conventional BCS case. Especially, for 
•q = 0.66, the uniform superconductivity does not appear 
at T = 0, but LOFF phase starts at T = 0. For suf- 
ficiently weak coupling, the uniform superconductivity 
disappears at any temperature, but the LOFF phase can 
survive at low temperature. In FigJ7]we show the neutral 
LOFF state for r\ — 0.63 where the uniform superconduc- 
tivity starts to disappear. The temperature behavior of 
the order parameter is similar to that in the conventional 
BCS case. 

We conclude that when LOFF phase is taken into ac- 
count, the phase diagram [411 - I45T ] of neutral quark matter 
is significantly changed. 
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FIG. 6: The pairing gap A in 2SC/g2SC phase (dashed line) 
and LOFF phase (solid line) and the LOFF momentum q as 
functions of T for two coupling values r\ — 0.70 and 0.66. 
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FIG. 7: The LOFF pairing gap A and LOFF momentum q 
as functions of T at extremely weak coupling r\ = 0.63 where 
the 2SC/g2SC phase starts to disappear. 



B. Gluonic Phase 

The gluonic phase is energetically favored at strong 
coupling where only the 4- 7th gluons suffer instability 
and may appear at intermediate coupling too(23. l40j|. In 
the gluonic phase, off-diagonal gluon condensate or the 
spontaneously generated off-diagonal Nambu-Goldstone 
current is nonzero. Since the 4- 7th gluons form a complex 
doublet, we can introduce only the 4-th gluon condensate 

P = g(A 4 ) = (Vip 4 ). (66) 

Including this condensate, the thermodynamic potential 
of the system becomes Q(A,p). Following the same pro- 
cedure used above, we can prove: 

(1) Thc neutral gluonic phase with p ^ exists only in 
the temperature region T < T 4 where the Meissner mass 
squared m| is negative. 

(2) The gluonic phase has lower free energy than the uni- 
form phase in the region T < T4. 

(3) At T — T4, the value of p approaches to zero continu- 
ously. 

The proof is quite similar to that for the LOFF phase. 
We need only the potential curvature definition of the 
Meissner masses squared [HJ 

d 2 n/dp 2 \ p=Q = m 2 4 /g 2 (67) 

and consider the fact that the gap equation for p = 
(0, 0, p) can be generally expressed as 

dn/dp = pK{p) = 0. (68) 

Recently, the numerical calculation on neutral gluonic 
phase is presented at T = in [40(. 



V. SUMMARY 



We have investigated the stability of neutral two fla- 
vor color superconducting quark matter at finite temper- 
ature. The main conclusions are: 

(1) There exists a temperature window below and close to 
the critical temperature of the superconductivity where 
the uniform 2SC/g2SC phase is stable, namely the Meiss- 
ner masses squared and the gap susceptibility are both 
positive. The Meissner mass squared is positive at tem- 
peratures T4 < T < T c for the 4-7th gluons and at 
T 8 <T <T C for the 8th gluon. 

(2) In a wide range of coupling in the diquark chan- 
nel, we have T4 > T§, the introduction of LOFF phase 
can not completely solve the problem of chromomag- 
netic instability at finite temperature, since in the region 
T s <T <T 4 only the 4- 7th gluons suffer instability. 

(3) Thc LOFF phase can exist only below the turning 
temperature T% where the Meissner mass squared for the 
8th gluon changes sign and the LOFF momentum ap- 
proaches to zero. Similarly, the gluonic phase can exist 
only below the turning temperature T4 where the Meiss- 
ner mass squared for the 4- 7th gluons changes sign and 
the off-diagonal gluon condensate vanishes. 

(4) When the LOFF phase is taken into account, the 
strange temperature behavior of the uniform color super- 
conducting order parameter Q disappears and the corre- 
sponding phase diagram of neutral quark matter (4lTj45j 
is significantly changed. This situation is quite like the 
recent studies for non-relativistic Fermi gas with popula- 
tion inbalance [32l - l34l | . 

For further investigation, one needs to make a detailed 
calculation for the phase with gluonic condensation or 
Nambu-Goldstone current and to check the chromomag- 
netic stability of the LOFF phase and gluonic phase. 
Since the LOFF momentum q and the gluonic condensate 
p approach to zero continuously at the turning tempera- 
tures T$ and T4, one can at least expand the thermody- 
namic potential in terms of q and p in the neighborhood 
of Ts and T4. We defer the research in this direction to 
be a future work. 
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Appendix A: Functions A, A' , B,C, D, H and J 

We list in this appendix the functions A, A' , B,C, D, H and J used to express the Meissner masses squared (|38|) 
and (|39j) for the 4-8th gluons and photon. Using the coefficients C± ± (p) and C± T (p) defined in [9( and taking the 
trick of replacement po — > — pa in calculating the Matsubara frequency summation, we can prove that only 7 of the 
coefficients are independent, 

A', C 12 _ = Cf_ = C 12 + = C 21 + = B, 
C 12 + = Cf + = C 12 _ = C 2 }_ = D, 

: J (Al) 
with the explicit expressions of A, A', B, C, D, H, J as functions of p = |p|, 

A(p) = u 2 _v 2 _ ±1^± 1 + u\vl tW> + l ( f' ) - 1 + uif(E 2 ) + vif( El ) + u%rm + </'(£ 4 ), 



rill 
°++ 




= A, 


r22 

°++ - 


C 22_ 




= c n + 


= cf_ 


= C 2 \ 


= c, 


r44 
°++ 


= C 44_ 




cf_ = 





A '(p) = utvt J -^ + u\v 2 + J ^' - +utf'(E 1 )+vtf(E 2 ) + u\fm + vlf'(E 3 ), 

^A A 

B{p) = u_v_ — — r (/ (J5i) + / {E 2 )) - u%v% (/ {E 3 ) + f (£4)) , 

E A E A 



C( P ) 
Dip) 
Hip) 



(E A ) 2 -E^E+ f(E 1 )+f(E 2 )~l {EX? - E^E+ f(E 3 ) + /(£ 4 ) - 1 1 
(^) 2 -(^) 2 ^ (^a) 2 - (EX) 2 E+ p 

A 2 + /(£ 3 ) + /(£ 4 ) - r 



(£ A )2-(£+) 2 V £ A 
ff(E 6 )-f(E 2 ) , f(E s )-f(E 1 )\ , 2 ff(E 6 )-f(-E 1 ) , /(£ 8 )-/(-£ 2 ) 



U 2 _ I — 1 + r 



2 



\ £e — £-2 -E-8 — Ei J \ Eq + £•! £/8 + £ 

2 f f(E 5 )-f(E 3 ) f(E T )-f(E 4 ) \ 2 ff(E 5 )-f(-E 4 ) , /(£ 7 )-/(-£ 3 ) 
u, ! + f. 



J(p) = u 2 _ 



V £5 — £3 Ef — Ei ) + \ £5 + £4 £7 + £3 

ff(E 5 )-f(-E 2 ) , f(E r )-f(-E 1 )\ , 2 f f(E 5 )-f(E 1 ) /(£ 7 )-/ (£ 2 ) 



\ -C/5 ~h -^7 "h -El / V ^5 — -^l -^7 — 



2 ff(E 6 )-f(-E 3 ) /(£ 8 )-/(-£ 4 )\ 2 ff(E 6 )-f(E 4 ) , f(E s )-f(E 3 )\ , 2 



+ V S 6 + £ 3 £ 8 + ^4 / + V ^6-^4 E 8 -E 3 J p : v ; 

where the coherent coefficients w 2 ^ and w j_ are defined as itj_ = (1 + E±/EX) /2 and u| = (l - E^/Ef) /2. Note 
that we have added the terms 1/p to C and 2/p to J to cancel the vacuum contribution. In this way the Meissner 
masses squared are guaranteed to be zero in the normal phase with A — 0. 
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